program squarewell
	use fs_FindSolution_module_q5b
	implicit none
	real, parameter	::	k = 26.2468512
	real			::	v
	real			::	a
	real 			:: 	root
	!for a = 0.3
	a = 0.3
	do v = 0.2, 20, 0.2
		fs_k = k*v*a*a
		root = fs_bijection(0.0001,fs_k**2/2.0,0.0001)
		print *, v, root/a
	end do
end program
